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ABSTRACT 

Among other things, studies of the formation and evolution of planetary systems currently draw on 
two important observational resources: the precise characterization available for planets that transit 
their parent stars and the frequency and nature of systems with multiple planets. Thus far, the study 
of transiting exoplanets has focused almost exclusively on systems with only one planet, except for 
considering the influence of additional planets on the transit light curve, mostly through transit timing 
variations (TTVs). This work considers systems where multiple planets are seen to transit the same 
star and concludes that such "multi-transiting" systems will be the most information-rich planetary 
systems besides our own solar system. Five new candidate multi-transiting systems from Kepler have 
been announced in Steffen et al. 2010, though these candidates have not yet been fully confirmed as 
planets. In anticipation of the likely confirmation of multi-transiting systems, we discuss the value of 
these systems in detail. For example, proper interpretation of transit timing variations is significantly 
improved in multi-transiting systems. The true mutual inclination, a valuable probe of planetary 
formation, can also be well determined in certain systems, especially through Rossiter-McLaughlin 
measurements of each planet. In addition, such systems may undergo predictable and observable 
mutual events, where one planet crosses over the other, which allow for unique constraints on various 
physical and orbital parameters, particularly the mutual inclination. 

Subject headings: stars: planetary systems, techniques: photometric, techniques: spectroscopic, 
eclipses, occultations 



1. INTRODUCTION 

Even as transiting extra-solar planets are the most 
information-rich planets orbit ing other stars (Charbon- 



neau et al. 2007 



Winn 20091, systems where multiple 



planets transit their parent star will be the most informa- 
tion rich planetary systems around main sequence stars 
after our own solar system. Though rarer than systems 
where one or no planets are seen to transit, these systems 
provide an exciting opportunity to learn more about the 
formation and evolution of planetary systems. 

As of May 2010 there were no published examples of 
systems with more than one planet known to transit. In 
the last year, two transiting planets were shown to have 
an additional companion well-characterized from radial 
velocity observations: the super-Earth CoRoT-7b (Legerl 
|et al. 2009; Qucloz ct al. 2009) and the hot Jupiter HAT^ 
P-13b (Bakos c t al.,2009^ . HAT-P-13c has a much larger 
orbit than planet b, and it is not known whether it tran- 
sits. CoRoT-7c is in a 3.6-day period orbit, which is 
not transiting. Both of these systems are indicative of 
the trend that hot Jupiters tend to be the only mas- 
sive object close to the star, while hot Neptunes and 
super-Earths show the opposite trend, tending to clus 



ter in multiple planetary systems (e.g., Lo Curto et al. 
2010| which h as been predicted theoretically ( |l'erquem 
& Papaloizou||2007 ). Though this inference is subject to 
observational bias, there are many examples from radial 
velocity surveys of close-in multi-Neptune systems with 
near circular orbits, which are quite likely to have low rel- 
ative inclinations as well: coplanar analogues of the 61 
Vir and HD 40307 systems that are inclined within 0.6° 
and 1.5° degrees of the line-of-sight, respectively, would 
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show all three pla nets transiting; HD 40307b is known 
not to transit (Gillon et al.||2010b[ ). 

The last decade of radial velocity observations has 
shown that the best way to find planets is to look for ad- 
ditional planets around stars with known planets. This 
applies to transiting planets as well: systems with tran- 
siting planets are much more likely t o contain additional 



transiting planets than random stars (Gillon et al. 2010a 



Beatty fc Seager|2010 1 . This is an extension of the early 



realization that if stellar binaries contained planets in the 
same plane, then eclipsing binaries are an ideal place to 
search for transiting planets (e.g., Borucki & Summers 
1984 Schneider fc Chevreton|1990 K Several groups have 



searched for additional transiting planets aro und systems 



with known pla nets, including MOS T (e.g., Croll et al 



[20071), EPOXI dBallard e"ral[2010|, and Super- WASP 



(Smith et al. 2009[ ). Unfortunately, these surveys have 



limited observations and statistical analyses show that 
they are seldom sensitive to planets with periods longer 
than ~10 days, implying that only very close additional 
companions would be detectable. 

Observations by Kepler have ushered in a new era of 
exoplanet studies, with the recent announcement of hun- 
dreds of cand i date t ransiting exoplanets in Borucki & the 



Kepler Team ( [2010 ), with an estimated false positive rate 
that implies that in the first 43 days of operations, Ke- 
pler discovered about as many new exoplanets than all 
previous methods combined over the last 20 years! Fur- 
thermore, this an nouncement was accompanied by |Stef-| 



fen et al. (2010 hereafter SIO) with describe in detail 
5 new systems where more than one planet candidate is 
seen to transit. Though these have not been confirmed as 
planets an d som e of these candidates could be false pos- 
itives (see S3.1), they demonstrate the ability of Kepler 
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to identify and characterize multi-transiting systems. As 
expected, the systems reported are mostly composed of 
smaller planets in long-period orbits. Most of the sys- 
tems appear to be near mean-motion resonances and one 
system has three planetary candidates. It is important 
to note that the 5 systems reported are only a subset 
of multi-transiting systems identified by Kepler to date 
(SIO). 

With their high sensitivi ty and long observ ing base- 



lines, the CoRoT mission (Baglin et al 



Kepler Space Telescope ( Borucki et al 



2006) and the 



0| are far more 



likely to find a system with two or more planets detected 
in transit around the same star. The lack of additional 
transiting planets around the reported CoRoT systems 
seems to indicate that multi-transiting systems are only 
<10% as common as tr ansiting planets, in accordance 
with previo us est imatesj Holman & Murray 2005 Fab- 
rycky |2009[ [Holman 20lo| . Without full disclosureo! 
the Kepler population, the observed frequency of Kepler 
multi-transiting systems is not clear, though they defi- 
nitely appear to cont ain more than ^1% of the detected 



planetary candidates (Borucki & the Kepler Team'2010l). 

In this work, our goal is to focus on understanding 
systems in which multiple planets are already known to 
transit and for which there will be significant signal-to- 
noise in the transit light curves of both planets, such as 
the systems presented by SIO. We focus on further an- 
ticipated results from the Kepler, which is currently re- 
turning ultra-prec ise light curves from over 150000 stars 
( |Koch et al.]|2010| . That Kepler would be likely to find 
and interpret su c h sys tems was anticipated as early as 
Koch fc Bor"^lAi| ( [19961 ). 

I'his paper is a systematic investigation of the value of 
multi-transiting systems. In S|2| we define our terminol- 
ogy, review the study of multiple planets around other 
stars, and consider the geometric probability of multi- 
transiting systems. Specific aspects of the value of multi- 
transiting systems in discussed in f|3l In f|4| we describe a 
full numerical-photometric model that can produce self- 
consistent light curves of multi-transiting systems. Ad- 
ditional methods for studying multi-transiting systems, 
including mutual events and Rossiter-McLaughlin, are 
discussed in We summarize with specific conclusions 
in ^ 

Much of the following discussion is relevant not only 
for systems with two planets in different orbits, but 
could also be directly or indirectly applied to systems 
with satellites (exomoons), trojans, transiting planets 
in eclipsing binary systems and "multi-eclipsing" stellar 
triples. The photometric effects of some of these have 
been modeled by other authors. 

2. BACKGROUND 

2.1. Terminology 

Systems where more than one planet is seen in tran- 
sit do not yet have a standard terminology, so we find 
it valuable to first define some terms. ^ Multi- planet sys- 
tems are stars orbited by more than one planet, whether 

^ We also note that there is no standard term for a triple star 
system where eclipses between more than one pair are occurring; 
eclipsing triples can refer to a regular eclipsing binary with a non- 
eclipsing third star. By extension with the terminology we use, 
these could be called multi-eclipsing triples. 



or not any are seen in transit. A multi-transiting sys- 
tem (MTS) is one in which more than one of these plan- 
ets (currently) transits the parent star and these are the 
subject of this paper. This can be made more specific, 
e.g. double-transiting system, triple-transiting system, 
etc. In most cases, not all of the planets will transit, 
so you can readily have a double-transiting triple system 
or a single-transiting quadruple system, etc. A mutual 
event refers to an overlap between two of the planets as 
seen from Earth in general and the terminology we use 
for these events (e.g., double transits, etc.) is described 



in J 5.1 



2.2. Transiting Planets with Companions 

Most of the literature dealing with multi-planet tran- 
siting systems is focused on a transiting planet perturbed 
by an unseen second planet causing non-linear varia- 
tions in the observed central transit time that cannot 
be explained by a purely periodic Keplerian orbit (e.g., 
Miralda-Escude[[20g2| |Agol et al.|[2005| [Holman fc Mur-| 
ray|2005[ ). Since transit times can with good observations 
be measured to within a few seconds^, in the best cases, 
the presence or lack of transit timing variations (TTVs) 
is a powerful constraint on any non-Keplerian effect, such 
as an additional planet. 

Though transit timing anomalies can be very large, as 
much as several minutes in the case of strongly resonant 
systems'^, no announced planets have been detected us- 
ing TTVs. The systems in SIO only cover 43 days of 
data and no significant TTVs are reported. A Monte 
Carlo analysis of expected TTVs over the course of the 
Kepler mission indicates that these systems, especially 
the near-resonant ones, could show strong TTVs in the 
future (SIO). We note that if these systems are n earl y 
coplanar, as suggested by geometric arguments (^2.3), 



the eccentricities estimated by SIO may also be low, as in- 
clinations and eccentricities are usually dynamically cou- 
pled, suggesting that the actual TTV signals will be on 
the low end of the expectations. 

At the current time, there are a few hints of possible 
varia bility due to unsee n companions in the literature 



(e.g., Maciejewski et al. 2010), but a strong clear TTV 



signal has yet to be demonstrated. Furthermore, the in- 
version process of taking transit times and inferring the 
properties of the perturbing planet appears is quite dif- 
ficult in general (Veras & Ford 2 010, submitted), though 
not imposs ible (Ncsvor ny fc Bea ugc 2010; N esvorny fc] 
Morbideni||2008^ .Meschiari fc Laughlin^ 2010). Besides 
a computationally difficult inversion prooTem, the TTV 
signal is often degenerate, with multiple possible config- 
urations that each fit the data, even i n systems where 
TTV s are measured quite accurat ely ( [Ford fc Holman 
20071 [Meschiari fc Laughlin|[2010[). This can be com 



plicated even further when allowing for more than one 
perturbing planet (SIO). 

Considering systems where both planets are seen to 
transit helps the inversion problem immensely by con- 

^ Measuring central transit times to within a few seconds is com- 
parable to measuring orbital phase position to within 1 arcsecond. 

^ Note that near-resonant systems also have large TTVs (Veras 
& Ford 2010, submitted) and that the apparent period of a reso- 
nant planet (during a small fraction of the libration cycle) can be 
different from the true period by a couple percent, so the ratio of 
the apparent periods is not the only distinguishing factor. 
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straining the period, phase, and size of the perturb- 
ing planet directly, without degeneracy. These param- 
eters are enough to get good estimates on the perturbing 
planet's semi-major axis and mass, the primary deter 



minants of the TTV signal (Agol et al. 



2005 



Holman 



& Murray 2005). Furthermore, as both planets perturb 



each other, there can be two sets of interdependent TTVs 
that must be self-consistent, which is particularly helpful 
in ruling out additional unseen perturbers (SIO). We dis- 
cuss the excellent value of such systems for constraining 
all the planetary system parameters in more detail below. 
Of course, using TTVs alone will still be needed in the 
case where the second planet does not transit and/or the 
second planet is small enough to be undetectable (ini- 
tially) in transit. This latter case might occur when the 
planets are in resonance, since this significantly amplifies 
the TTV signal. 

There is also some information in the changes in tran- 
sit durations or sh apes due to additional planets, but 
( Ragozzine fc Wolf||2009| ) conclude that, in the majority 
ot expected exoplanet systems, the signal due to tran- 
sit timing anomalies is usually stronger than that due 
to variations in the transit shape, unless the planets are 
in specific non-resonant orbits with large mutual inclina- 
tions. We confirm this notion and have found, for ex- 
ample, that the amplitude of transit duration variations 
(TDVs) for HAT-P-13b as induced by HAT-P-13c are 
less than half the size of the transit timing anomalies, at 
any relative inclination. Note also that the accuracy of 
transit duration measurements is usually less than half 



as go od as transit timing measurements (Carter et al. 



20081 and that transit durations are much more sensi- 



tive to possible systematic effects such as observing in 
different filters. Although "transit shaping" is harder to 
detect than transit timing, when it is reliably detected, 
it can have useful value for determining the properties 
of the planetary system. Note that when the star ex- 
periences significant acceleration during the course of a 
transit, TTVs and TDVs may not be entirely well defined 
(see Sji] below). 

Beyond SIO, a few studies have considered systems 
where more than one planet is found to transit. Anoma- 
lies in photometry from TrES-lb from Hubble and the 
ground were co njectured to b e possibly due to an ad- 



ditional planet (Rabus et al. 2009| ), although the au 



thors t hemselves acknow l edge this as an unlikely pos 



Dittmann et al. (2009) argue that such signals 
1 better explained by crossing starspots, espe- 



sibility. 
are muc 

cially for TrES-1, which is known to be photometrically 
active. Authors have noted that nearly coplanar multi- 
planet systems are relat ively likely (>10%) to have mul- 
tiple planets tr ansiting fSchneider"2004! l E^brycky||2009 
Holman & Murray 2005). Schneider (200?) also quanti- 
tatively showed that confusion between which planet is 
tr ansiting at any particular time is extremely unlikely. 



Minniti (2005) submitted a manuscript to arXiv that 
was never published, with a claim of observing a transit 
of a second resona nt planet in OGLE -TR-111 system; 
new observations by Adams et al. ( 2010 1 have shown that 
there are no significant TTVs in this system. 

Though SIO report five new candidate multi-transiting 
systems, as far as we are aware, there is no known in- 
stance of an object that transits or eclipses one or both 
stars in an eclipsing binary (D. Fabrycky, G. Torres, pers. 



comm.), though many eclipsing binaries are known to 
have additional companions (i.e., from the stellar equiv- 
alent of TTVs) . Eclipses of the same star by two stellar 
companions are a priori unlikely since the large masses 
of three stars imply that gravitational stability is only 
possible if the triple is hierarchical, requiring the outer 
star to have a large orbit and making eclipses observable 
from Earth relativ ely improbable. Despite early searches 
(e.g., Deeg et aL][r9 98), to date, no transiting pla nets 
have been touna ar ound eclipsing binaries, though |A1-| 
menara et al. ( 2009[ ) propose an algorithm for detecting 
such systems. 



2.3. Geometric Probability of Observing 
Multi- Transiting Systems 

From radial velocity observations, there is growing ev- 
idence that hot Neptune/Super-Earth planets with short 
periods (P < 50 day s) are common, wi th a reported fre- 
quency around 30% ( Mayor et al.|2009 |. It is also appar- 
ent that hot Neptunes tend to cluster in multiple systems 
(Lo Curto et al. 2010) and that these could potentially 
show multiple planets transiting (SIO). This indicates 
that the probability of finding multi-transiting systems, 
especially using with the precise photometry and long 
baseline of Kepler^ will be relatively high. Can we place 
a quantitative estimate on the number of multi-transiting 
systems we expect? 

While a full answer requires a detailed model of planet 
formation and evolution (e.g. 



Ida & Lin 2010), we can 



straightforwardly address the geometric aspect of this 
question. In particular, we consider the following. Sup- 
pose that all stars that have a transiting planet (1) also 
have an additional planet (2) at a true mutual inclina- 
tion of cos(/)i2 = cosii cosi2 -I- sinii sini2C0s(ni — 0.2)- 
(These angles are defined in ^j4j) What is the probability 
that, given planet 1 is transiting, that planet 2 will also 
be transiting, as a function of (j) and the semi-major axes 
of the two planets? 

To answe r this geometric question, we turn to the 
method of (Borucki & Summers 1984), later repeated 
by other authors, who consider the traction of the ce- 
lestial sphere centered on the parent star where distant 
observers would see the planet (s) transit. A model of 
this is shown in Figure [T] There is a certain fraction 
of the celestial sphere where distant observers would see 
no planets transit, the inner planet transiting, the outer 
planet transiting, or both planets transiting. This de- 
pends both on the semi-major axis (in stellar units) of 
each planet, as well as the true mutual inclination. There 
are two cases: 1) the mutual inclination is low, so that 
the entire region that can see transits by the outer planet 
can also always see transits by the inner planet and 2) the 
mutual inclination is high enough that observers must be 
near the line of nodes (the intersection of the two orbital 
planes) to observe both planets transiting. 

In the low mutual inclination case, the probability that 
both planets are seen to transit is just i?*/a2, the transit 
probability of the outer planet. In the high mutual incli- 
nation case, an estimate of the probability can be made 
by ignoring the spherical geometry and calculating the 
area of the parallelogram on the celestial sphere near the 
line of nodes where a distant observer would see both 
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Fig. 1. — Illustration for determining frequency of multi- 
transiting systems. The celestial sphere, centered on the parent 
star, can be divided into regions where no planets are observed 
to transit (white) , only one planet transits (lighter shade of gray) , 
and both planets transit (darkest shade of gray). (Only parts of the 
closer hemisphere of celestial sphere are shaded in this illustration.) 
In a latitude band centered on the planet's (current) orbital plane, 
distant observers would see the planet transit; the size of this band 
is given by ^ , in the case of circular orbits, with the more distant 
planet having the narrower band. An additional planet is inclined 
at the true mutual inclination </>. In the overlap region (dark gray), 
both planets are seen to transit. The apparent on-the-sky inclina- 
tions for the planets depend on the exact location of the observer 
in this region. The area of this region, divided by the surface area 
of the sphere, gives the probability of seeing both planets transit. 
See discussion in text and Figure [5] 



planets transit (Figure [ij, which gives 



p(l&2) 



0x02 Sine 



(1) 



where p is used for probabihty and "1&2" means both 
planets transit. Both geometrically and from continuity, 
it can be seen that the inclination which divides these 
regimes is (/icrit — sin</)crit — although this is ap- 
proximate as it also depends somewhat on the spherical 
geometry of the system. 

These are the probabilities of observing both planets 
in transit given no a priori information. If we know that 
one planet is transiting, then we can calculate the proba- 
bility of the other planet transiting given that one planet 
is known to transit. That is, we can now answer the pro- 
posed question, what is the probability that the planet 2 
will transit if planet 1 is known to transit. The approxi- 
mate answer is (see also Figure [2]): 



p(l&2|l 




if «! > 02 and 



if, 



ai 



a-[ < a-\ and 6 ^ — 
— ^ if (/) > ^ 

02 Sin ^ ^ ai 



These analytic expressions are shown with the results 
of a simple Monte Carlo model in Figure [2] In this model 
two planets are set up with mutual angles 0, and are ob- 
served from 10*^ lines of sight placed randomly on the 
celestial sphere. The model then calculates the appar- 
ent on-the-sky inclinations, and then gives the probabil- 
ity that both planets are seen to transit if the planet 
1 is seen to transit. We ignore the size of the planet 



and non-circular orbits; these can be incorporated by 
adjusting the size of the star or, for eccentric orbits, us- 
ing the "effective" semi-major axi s that gives the same 
star-planet distance r = aj ( |Winn|2010[ ). We also 



assume the planets to be on fixed orbits, though interac 
tions over many orbital periods can result in precessing 
orbits that go in and out of transit (e.g., Schneider & 



Doyle 1995 Shankland et al. 2006). The two regimes 



can be seen and the overplotted analytical expressions 
provide a reasonable approximation; differences between 
the analytical and numerical calculations result from not 
strictly incorporating the spherical geometry and dimin- 
ish with increasing semi-major axes. Our result is com- 
parable to, but more general th an, the derivation of the 
double-tra nsiting probab i hty b y |Koch fc Borucki ' 
and Holman fc Murray] (|2005[) b 



19961 
tEe 



jecause we inclu' 
appropriate random distribution of the inclination of the 
transiting planet. These analytical expressions suggest 
that the probability of seeing an both Venus and Earth 
transit, given an exact analogue, is ~0.06%, though if 
the Venus analogue is known to transit, then the proba- 
bility of the Earth analogue also transiting is 8.4%. For 
the Super-Earth system triple HD 40307, if the mutual 
inclination is low (;$1.5°), then the probability that all 
three transit is essentially just the probability that the 
outer planet transits, ~2.6%. The probability of observ- 
ing systems like those recently announced is given in SIO. 

We have verified that the above results are consistent 
with the random 3-dimensional placement of the Laplace 
plane and a random orientation of the two planets with 
respect to that plane. Assuming a two planet system, 
the inclination distribution (from the plane of the sky) 
of the second planet can be determined by: 



cos i2 — cos ii cos (j) + sin i i sin </> cos A 



(2) 



where A is a randomly chosen angle (and inclinations 
are in the observer's reference frame, i = 90° is edge- 
on). The distribution of 12 — ii is not Gaussian. In- 
deed, it is peaked at the extremes: «2 is most likely to be 
ii + 4> OT ii — (f), though there is still significant probabil- 
ity throughout the range. (Inclinations, 12, greater than 
90 degrees are equivalent to 180° — 12 ■) Note that this 
accounts for the full three-dimensionality of the problem, 
rather than assuming a typical on-the-sky inclination dif- 
ference (which is only weakly related to the true mu- 



tual inclination) with some distribution (Beatty & Gaudi 
2008 1 . The above discussion assumes there is no system- 
atic relationship between the orbits of the planets, but 
this is not necessarily true when the mutual inclination 
is in the Kozai regime (40° ^ </> ^ 140°), where the secu- 
lar Ko zai effect can affect the or ientation of the planetary 
orbits ( |Kozai|1962] |Lidov|1962| . In particular, the Kozai 
effect acts to bring the periapse of inner planet to occur 
when it is maximally out of the plane of the outer planet; 
this would effect the probability of observing this system 
as double-transiting and should be kept in mind. 

Based on these results we can draw the following infer- 
ences. The fraction of double-transiting systems remains 
relatively high (^) even for non-zero mutual inclinations 

(0 < ^). The true number of multiple systems with one 
planet transiting is usually a few to several times larger 
than the number of systems where both planets are seen 
to transit, even if all systems are nearly coplanar. This 
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Fig. 2. — Based on the geometric model presented in Figure^ 
and described in the text, the probabihty that a two-planet system 
would show two planets transiting, assuming that the first planet 
transits, can be calculated as a function of the true mutual incli- 
nation (which is only weakly related to the apparent inclination 
measured by transit durations) ranging from coplanar (0°) to per- 
pendicular (90°); the curves are symmetric around <p = 90°. It 
is important to note that the calculated probability assumes that 
all systems are two-planet systems, i.e., is a geometric argument 
only and does not say anything directly about the astrophysical 
likelihood of double-transiting systems. The solid lines give the 
probability (denoted 1&2 — 1) for various systems as determined 
by a Monte Carlo model, which accounts for a random viewing an- 
gle. The planetary systems in question are identified by the legend: 
the first number is the orbital period in days of the planet that is 
assumed to transit and the second number is the orbital period 
in days of the planet whose transit probability is being assessed. 
For example, a system with a 3-day transiting planet and an ad- 
ditional planet at 9 days is shown with asterisks. An analytical 
estimate, divided into two inclination regimes, is described in the 
text and given by the dashed line for the system with planets at 
3 and 9 days. We have also shown the probability that a 1-day 
planet transits assuming that a 10-day outer planet is known to 
transit (diamonds); at low mutual inclinations, transits of the 1- 
day planet are assured. Note that the reciprocal of the probability 
given above is the number of similar systems that would be ex- 
pected in the full population where the second is present but not 
seen to transit. 



reinforces the fact that the search for TTVs is still an im- 
portant method for probing multiple planetary systems. 
In particular, the types of systems where multiple plan- 
ets are seen to transit will be strongly biased against 
those with large (> 5°) mutual inclinations. Even if 
the fraction of double to single systems detected by Ke- 
pler is comparable to the values expected from nearly 
coplanar orbits, this does not necessarily imply that all 
systems have low mutual inclinations. One should com- 
pare the distribution of other orbital parameters (e.g. 
eccentricities, spin-orbit alignments) of multi-transiting 
systems with, for example, radial velocity multiple sys- 
tems as geometric biases imply that multi-transiting sys- 
tems are drawn from the low-inclination subset of multi- 
ple planet sy stems. Statistical t e chniqu es, such as those 
employed in Fabrycky & Winn (2009), can make more 
rigorous statements about the distribution of mutual in- 
clinations with strong implications for planet formation 
as discussed in §3.3| 

In systems with more than two planets, secular evolu- 
tion of the mutual inclination is possible and n-body in- 
teractions (i.e., as characterized by second-order Laplace- 
Lagrange theory) are relevant to the evolution of the 
planetary inclinations and nodal longitudes. This im- 



plies that the orbital plane may be constrained in ways 
not considered here. 

3. THE VALUE OF MULTI-TRANSITING SYSTEMS 

Multi-transiting systems are valuable for understand- 
ing the physical and orbital parameters of distant plan- 
etary systems for at least two reasons. First, since both 
planets transit, we have a wealth of information about 
each planet that is not known for non-transiting plan- 
ets, such as the radius and inclination, or is much more 
precisely determined, such as orbital period and phase. 
This already makes multi-transiting systems quite valu- 
able. The second valuable aspect is that both planets 
are known to transit the same star and these planets are 
often close enough to affect one another's orbit. This syn- 
ergy enhances the value of these systems and the follow- 
ing discussion attempts to focus on this aspect of multi- 
transiting systems. In particular, we consider the value 
of such systems in eliminating false positives, measur- 
ing relative and absolute masses and radii, determining 
mutual inclinations, providing additional dynamical in- 
sights, and for use in comparative planetology. 

3.1. Eliminating False Positives 

Many transiting systems, especially those discovered 
by Kepler, will be around faint, active, and/or rapidly 
rotating stars for which it is difficult to get abundant and 
accurate radial velocities. Observation time to confirm 
the hundreds of good planetary ca ndidates from Kepler 
( [Borucki fc the Kepler Team|2010"| will also be very lim- 
ited, esepcialiy given the observed preference for lower 
mass planets. In the case of single transiting planets, ra- 
dial velocity confirmation is generally used to show that 
the light curve is not created by any other astronomical 
phenomenon and to measure the minimum mass of the 
putative companion. Astrophysical signals that mimic 
the light curve of transiting planets (i.e., producing a 
~1% drop in light on several hour timescales) have been 
a serious issue for ground-based transit surveys. These 
astrophysical "false positives" are usually very common 
and include, for example, an eclipsing binary contami- 
nated by the third light of a brighter (possibly related) 
star, a dwarf star orbiting a giant, or a grazing eclipsing 
binary"* . With its astounding photometric and astromet- 
ric precision, Kepler can eliminate many, but not all, of 
these false positives even before spectroscopic investiga- 
tion (Batalha et al. 2010). In any case, the hypothesis 
that the light curve is caused by transiting planets can 
either be verified through radial velocity measurements 
that measure the mass or by eliminating all plausible as- 
trophysical false positives. This section will be concerned 
with the latter technique, since many transiting planets 
will not be able to be fully characterized with radial ve- 
locities due to astrophysical or logistical reasons. 

Multi-transiting systems are more difficult to mimic 
astrophysically than single-transiting systems, though it 
is not impossible. Some examples include, two eclipsing 
binary systems, either both grazing or contaminated by 
each other's light (or the "fifth" light of ano ther star). 
Such systems a re rare, but not unheard of (Pilecki & 
Szczygiel 2007). Multi-transiting systems can also be 



We will not discuss here false positives due to instrumental or 
systematic effects, which also need to be carefully ruled out. 
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comprised of a combination of false positives, planets 
around different stars, and planets around the same star. 
For this reason, each candidate transit signal should be 
investigated separately. Note that there are two steps 
in eliminating false positives in multi-transiting systems: 
showing that the candidates transit the same star and 
then investigating the false positives in the case of mul- 
tiple objects orbiting the same star. 

With precise photometry and astrometry, as is avail- 
able from Kepler^ several signs of false positives can be 
investigated: a secondary eclipse that is too deep to be 
caused by a planet; ellipsoidal oscillations indicating high 
masses; "Doppler photometry" or relativistic be aming 



the same star is the observation of TTVs, TDVs, or mu- 



from the stellar motion, that indicates high masses ( Loeb 
fc Gaudi||2003 ) ; multiple distinct stellar rotatio n or pul 



sation perio ds; variability typical of giant stars (Jenkins 
et al.|[20T0 ); astrometric centroid motion during eclipse, 



indic atmg a background eclipsing binary (Batalha et al. 
20101; and other methods. Most of these techniques ap 



pear to be used regularly by the Kepler team (SIO) and 
should be applied to each individual planet candidate 
with the requirement of consistent results, e.g., appar- 
ent astrometric motion due to crowding should go in the 
same direction for both candidates. We note that Kepler 
here does have one disadvantage compared to CoRoT in 
that it measures only mono- chromatic light curves and 
observations in multiple colors are an excellent way of 



ruling out false positives (Rosenblatt 1971 Leger et al. 



2009| . Observing even one transit ot' each candidate in 
a different filter is an excellent way to rule out false pos- 
itives, besides the added benefit of helping to break de - 
generacies involving limb darkening ( Colon & Ford 2009 ). 
Observations of the Rossiter-McLaughlin ett'ect can also 
be used to rule out false positives. 

There are a few techniques for eliminating false posi- 
tives that are possible only in multi-transiting systems. 
To show that the two planet candidates orbit the same 
star, one can check that the stellar properties inferred 
from the transits are consistent between each other and 
with the expectations for the stellar properties (e.g., 
based on spectral type). For example, the durations of 
the transits compared with the orbi tal periods can be 



used to estim ate the stellar density (Seager & Mallen- 
Ornelas 2003). This technique is not as powerful as might 



be expected because, though the stellar density can be 
measured precisely by a transiting planet on a circular 
orbit (if the impact parameter is well determined), al- 
lowing for a non-zero eccentricity introduces significant 
uncertainty in the stellar density measurement, which, 
to first order in eccentricity, is proportional to (1 -I- 3e) 
(SIO). There is a maximum eccentricity for each planet, 
assuming they are in orbit around the same star, from 
stability constraints. Therefore, if the inferred stellar 
densities from each planet are more different than allowed 
given the eccentricity constraints, it can be assumed that 
the candidates do not orbit the same star. Similarly, in- 
ferred limb darkening coefficients, if known accurately 
and without degeneracy, must be similar if the candi- 
date planets are transiting the same star. While these 
are important ways of identifying false positives, partic- 
ularly for high signal-to-noise planets, in practice, these 
properties alone do not seem sufficient for proving that 
candidates are indeed planets orbiting the same star. 
A much stronger indication that the candidates orbit 



tual events (^5.1). This is true for single-transiting sys- 
tems as well, but in multi-transiting systems if the TTVs 
are strongly consistent with perturbations from the other 
known transiting object, then this is a strong indication 
that both candidates orbit the same star. Care needs to 
be taken here to note that TTV measurements can be 
quite degenerate and/or affected by additional unseen 
perturbing planets. Interdependent TTVs that are con- 
sistent with each object perturbing one another provide 
clear evidence that both objects are in orbit around the 
same star. Whether TTVs are observed for one planet or 
more, if the masses of the planets may also be measured 
or constrained, this is equivalent to radial velocity mass 
meas urements and proves the planetary hypothesis (see 
^3.2 below). Indeed, once it is known, through whatever 
means, that both candidates orbit the same star, then 
a measurement or upper limits on TTVs will usually be 
sufficient to rule out the possibility that one or more of 
the candidates is stellar in nature. Orbital stability can 
also be invoked and will often be enough to rule out stel- 
lar masses. This is because stars are ~10^~^ times more 
massive than planets in general; brown dwarfs will be 
more difficult to rule out. Orbital stability can be used 
to identify false positives as well: if two candidates have 
very similar periods, then their orbits may be too close to 
one another to be stable, even with low planetary masses. 
These can be a chance superposition. 

It is important to note here that numerical analyses of 
multi-transiting systems, i.e. using the model described 
in f|4] are a valid and valuable method of "follow-up" to 
confirm the true planetary nature of these systems. 

Continuing along this theme, there are some expec- 
tations for the orbital spacing in multi-planet systems 
both observationally from the known RV planets and 
theoretically from planet formation simulations (though 
if the field of exoplanets has taught us anything, it is that 
planet formation theories may be incomplete). For ex- 
ample, planets in the same system are much more likely 
to lie near mean motion resonances through various for- 
mation and evolution effects, while planets around dif- 
ferent stars will only a ppear to be nea rly resonant by 
coincidence (compare to Goldreich|1965 ). Near- resonant 
architectures are already seen in the candidate multi- 
transiting systems of SIO and are a good indication that 
these systems are planetary in nature. System archi- 
tecture is especially convincing in multi-transiting sys- 
tems with more than two planets. Imagine a quadruple- 
transiting system that passes all the basic tests for false 
positives for each of the candidate planets. Any non- 
planetary scenario contrived to match the light curve 
will have ridiculously low probabilities: it is not reason- 
able to consider 8 equal brightness stars each in perfect 
eclipsing binaries with randomly selected periods having 
typical planetary spacing. You can then say with con- 
fidence that at least some of these candidates must be 
around the same star. Similarly, systems with apparent 
orbital periods that are so close together that stability is 
difficult may be an indication of a chance superposition 
of candidates around different stars. Also, given our un- 
derstanding on the mass dependence of multiple systems 
as discussed above, a system with multiple Neptune-size 
planets seems more likely than a system with multiple 
Jupiter-like planets. 
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Even after the planetary nature of a candidate is con- 
firmed, blending can be an issue for interpretation. For 
example, a transiting Jupiter blended with a brighter 
star will have the depth of a much smaller planet. How- 
ever, the full transit shape of a blended Jupiter and a 
true Neptune are not exactly the same (e.g., different 
ingress/egress durations) and Kepler is often powerful 
enough to resolve this degeneracy to some extent. 

Finally, we should note that systems which contain a 
false positive can still be scientifically interesting (e.g., a 
planet in an eclipsing binary). 

3.2. Dynamical Masses and Absolute Radii 

Kepler's third law relates the measured semi-major 
axis and period to the sum of the masses, but does not 
provide any information on each mass individually. In 
gravitationally interacting three-body systems, the de- 
generacies between the masses are broken and each of 
the mass ratios between the objects can be obtained in- 
dependently. This is valuable for several reasons. 

First, small mass ratios imply that the companions 
must be planetary in nature. That is, measuring mass 
ratios from TTVs is an extremely valuable method for 
rejecting false positives, especially when interdependent 
TTVs are simultaneously m easured, as discussed above. 

Second, Agol et al. ( 2005 ) point out that adding radial 
velocities to the transit timing variations gives dynamical 
masses for all of the objects in the system. Photometry 
alone combined with Kepler's third law measures the sys- 
tem parameters well, but the solution is invariant under 
transformations that conserve the stellar density, mass 
ratios, and radius ratios. In single-transiting systems, 
the mass ratio between the planet and the star must 
be provided by stellar RV measurements and the stellar 
mass is inferred through stellar modeling; often the the 
precision in absolute planetary parameters is limited by 
uncertainties in the stellar models. However, when TTVs 
accurately measure the mass ratio, the RV signal can 
be directly conve rted in the dynamical mass of the star 
( Agol et al.|2005 1. From the ph otometry, the stellar den- 



sity c an be precisely measured (Seager & Mallen-Ornelas 
2003 1 , so the radius of the star can also be determined 



m an absolute sense. Combining this with the precise 
radius ratios leads to very well determined radii for each 
object. When TTVs and RV are used in concert, there 
is no need for stellar models; indeed, such systems could 
be used as standards to improve stellar and planetary 
models. That the combination of TTVs and RVs yields 
dynamical m asses was confirmed numerica lly with simu- 
lated data by Meschiari & Laughlin (20101. 

It is also possible to solve for dynamical inasses without 
using radial velocities at all. If the light-travel timing off- 
set has a detectable TTV signal, it sets an absolute scale 
and measures dynamical masses without RVs. Generally, 
these shifts will be small unless the planets are massive, 
eccentric, and/or have strong TTV signals (Agol et al. 
2005; Veras & Ford 2010, submitted). Single-transiting 
planets have a small light-travel time effect during the 
course of the transit (Locb 2005) that is not detectable 
in pr actice even with Kepler photometry (Ragozzinc fc] 
Wolf 2009). However, multi-transiting systems provide 
both a precise clock (the inner transiting planet) and 
a well-known perturber (the outer transiting planet), so 
these are the systems that are most likely to have dynam- 



ical masses measurable from light-travel time TTVs; in 
these systems, dynamical masses are measured by pho- 
tometry alone. Note that the light-travel time effect has 
to be measured with high accuracy in order to provide 
a useful limit on the overall scale. It is also important 
to keep in mind the effects of potential unseen planets, 
particularly distant planets, which can confuse the inter- 
pretation of the light-travel time offset. Each system will 
need to be considered independently to see whether RV 
observations can easily provide a tighter constraint. 

The stellar radial velocity also has a small photo- 
metric effect due to relativistic beaming (at quite non- 
relativistic speeds) , which is known as photometric 
Doppler boosting ( Loeb fc Gaudi|[2003 1 . The photomet- 
ric amplitude of this signal is usually a few times the 
radial velocity signal divided by the speed of light, e.g., 
a large RV signal for a very hot Jupiter of ^1 km/s would 
result in a photometric amplitude of ~10 parts per mil- 
lion (ppm). On a faint Kepler star {V « 14, <t « 1000 
ppm/min), the photometric variation of such a best case 
system could be detected at the 3-a level in 21 days, but 
only in the absence of astrophysical noise. With planets 
that are smaller, longer period, and in the presence of 
astrophysical noise sources with comparable amplitudes 
(spotted stellar rotation, reflected light curves, ellipsoidal 
oscillations, etc.), this detection will be much more dif- 
ficult. Even so, it is worth noting that the combination 
of this effect with TTVs again could, in theory, provide 
a measurement of dynamical masses and absolute radii 
from photometry alone. 

3.3. Mutual Inclination 

The extent to which extra-solar planets are mutually 
inclined is a strong constraint on theories of planet for- 
mation and in understanding the evolution of plane- 
tary systems, whether by dis k-induced migration (e.g., 
iTerquem & Papaloizoul 2007 ), planet-plane t scattering 
{Chattcrjec et al. 2008/|Juri6 fc Tremain el [2008 ; Mat 

[Fabrycky fc 



sumura et al. 2iriTJp , Kozai ettects (e.g 
|Tremaiiic^2007j ) , or other mechanisms. Each of these will 
produce a different inclination distribution, and, when 
coupled with the eccentricity distribution and measure- 
ment of the Rossiter-McLaughlin effect (j |5.2[ ), promises 
the answer to many o f our questions about the origins 
of planetary systems ( Fabrycky 2009 ) . It is important 
to note here that observations of multi-transiting sys- 
tems do not immediately reveal the actual mutual in- 
clination of t he planets, whi ch is the quantity of great- 
est interest (!Fabrycky"2009). Even planets that have 
large mutual inclinations can lead to multi-transiting sys- 
tems in the correct orientation, though the probability 
of observing s yste ms shrinks rapidly with large mutual 
inclinations (S2.3). This implies that the number and 



distribution of multi=transiting systems, as an ensem- 
ble, will be extremely valuable for increasing our under- 
standing of the formation and evolution of planetary sys- 
tems (SIO). The importance of multi-transiting systems 
in constrai ning the mutual i nclinations was recognized 



early, (e.g., Rosenblatt 1971 Borucki & Summers 



1984 



Koch fc Borucki 1996^ 

I'he astrophysically valuable angle is the true mutual 
inclination cos(/>i2 = cos ji cos 12 -I- sin ii sinz2 cos(f2i — 



1)2) is a powerful indicator of past evolution, but it is 
not measured directly even in multi-transiting systems 
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(see ^ for definition of angles). It can only be weakly 
constrained as ii + i2 < </) < 180° — ii —12, with most 
of the likelihood at the extremes of this distribution, but 
significant possibility anywhere in that range. However, 
there are several techniques that can improve our knowl- 
edge of the actual mutual inclination. 

One way of measuring <f> is to resolve the or bit on the 
sky through astrometry (e.g., |McArthur et al. 



torowicz 



2010 



terferometry (jUnwin et al. 2008), p olarimetry (e.g . 



2009[), or direct imaging (Marois et al 



ik- 



I'his can measure the nodal angles directly and, in com 
bination with photometry, could yield good constraints 
on the mutual inclination. With current observational 
technology, these techniques are only possible in nearby 
and/or bright systems, and measurement of precise mu- 
tual inclinations is rare, though future instruments may 
make these observations possible. The mutual inclina- 
tion for the resonant planets in GJ876 was measured 
from an ex tensive and long set of radial velocity mea - 
surements (Correia et al. 2010 Bean & Seifahrt 2009), 
but this method is not generally applicable to other sys- 
tems and, although the true masses of the planets are 
now known, their radii and densities are not known since 
they are not transiting. 

Generally a measurement of absolute orientation on 
the sky is not possible, and so a more "indirect" measure 
of the mutual inclination will be required. For multi- 
transiting systems, since (p only depends on the mutual 
nodal angle, any measurement that can determine the 
relative nodal position (fii — ^,2) also measures the true 
mutual inclination; we will make use of this fact in our 
discussion below. It is already known that mutually in- 
clined systems undergo nodal precession which changes 
the observed impact parameters of the systems and sig - 
nificantly changes transit shapes ( Miralda-Escude|2002 ) . 
Mutual inclinations also can be constrained by '1"1' V s 



(Bakos et al. 2009 



which is most readi 



Nesvorny||2009[ [Payne et al.||2010 l 
y applied in multi-transiting systems 



With the knowledge of planetary M sin i from radial ve- 
locities, the lack of TTVs or TDVs and/or a require- 
ment on orbital sta bility can give an upper limit on the 
mutual inclination (Schneider 2004). We present new 



promising methods of measuring mutual inclinations of 
multi-transiting systems in ^ 

3.4. Comparative Planetology 

Anot her valu e of mu lti-transiting systems was pointed 
out by [Mimiiti ( 2005 ) : uncertainties in the stellar mass 
and radius propagate equally to all the planets, meaning 
that the inferred mass and radii of the planets (whose 
uncertainties can be dominated by uncertainties in the 
stellar properties) can be compared more directly. It is 
also exciting to note that planets in the same system 
are likely to have some comparable properties (e.g., age 
and perhaps metallicity) , so that comparison of internal 
models for planets in the same system could be a more 
powerful discriminant of interior processes than com par- 
isons between planets around different stars Guillot (e.g., 
2010). The effects of irradiation or tidal heating can be 



seen more strongly in contrast since many of the other 
variables are probably consistent from planet to planet 
in the same system. Furthermore, correlations between 
the planetary system architecture (where planets are rel- 
ative to one another) and the bulk composition of these 



planets as inferred from their densities will also be quite 
revealing, especially in systems with three or more plan- 
ets (SIO). In our own solar system, the rocky planets 
are interior to the gas giants which are interior to the ice 
giants, but exoplanetary systems (that have likely under- 
gone migration to increase their transit probability) may 
show interesting differences. 



From geometric constraints (^2.3), multi-transiting 
systems are likely to have relatively small orbital peri- 
ods, so that tidal interactions with the parent star could 
be non-negligible. Detailed orbital knowledge derivable 
from multi-transiting systems will help us understand the 
interplay between tidal forces and secular planetary per- 
turbations and potentially measure the Love number k2p 
(which probes the planetary interior structur e) and tidal 
dissipation parameter Qp for the inner planet ( Wu fc Gol 



dreich|2002[[Mardling|2007||Batygin et al.|200!j||Mardling 



2010) 



The Love number is a unique and novel probe of the 
planetary interior that provides additional information 
beyond th e planetary density as a measure of central con - 
densation (Ragozzine fc Wol f|2009| [Henning et al |2009 ). 
The aforementioned method tor measuri ng fc2p is indirect 
and only works for small inclinations (Mardling 2010) 
and with planets whose orbits are very well character- 
ized; multi-transiting systems are the best opportunity 
for these inferences. 

The tidal love number k2p can also be measured di- 
rectly by observing apsidal motion, a well-known prac- 
tice in the stellar binary community that has recently 
been applied to exoplanetary systems (Ragozzine & Wolf 
2009 ) . Since multi-transiting systems may be well char- 



acterized with dynamical masses not resulting from the 
direct application of Kepler's Law, it is plausible that 
some multi-transiting systems can be used to constrain or 
measure additional non-Keplerian parameters ^ including 



general relati vity or k2p for the interior planet (Ragozzine 
fc Wolf| 2009[). The unique determination of limb darken- 



ing by multiple planets in transit may also improve the 
possibility of measuring planetary asphericity, another 
probe of k2p, which is usually degenerate with limb dark 



ening parameters in single-tr ansiting syst ems (Barnes fc 
iFortney 2003; Rago zzine C Wolf 2009; C arter fc VVnm 
|2010| ). In all, we expect that multi-transiting systems 
will be very useful for studying diverse and unknown ex- 
oplanetary interior structure. 

We also note that, in systems where one planet is al- 
ready known to transit, single events from very long pe- 
riod planets (P ^ 1 year) can be interpreted accurately. 
These same events around random stars are less robust 
and less likely, thus multi-transiting systems provide a 
unique opportunity to probe the (very) rare distant plan- 
ets that transit their parent star. 

3.5. Additional Dynamics 

In systems where more than one transiting planet is 
known to be orbiting the same star, the durations of cen- 
tral transits by planets with circular orbits have a fixed 
ratio. Indeed, it is this fact that allows for a precise de- 
terminatio n of the stellar density from transit light curve 
modeling (Seager fc Mallen-Ornela s 2003). While allow- 

^ The best multi-transiting systems ma y also potentially serve a s 
tests of modified theories of gravity (e.g., [lorio fc Ruggiero|2009^ . 
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ing for eccentric orbits and non-central transits makes 
this more difficult to use as a technique for identifying 
false positives (SIO), once multiple planets are known 
to be orbiting the same star, the transit durations can 
be used the other way as a constraint on eccentricities. 
This is because the period ratio is well known, giving the 
semi-major axis ratio and hence the orbital velocity ra- 
tio. As this velocity is measured by transit duration, as 
long as the transits have high signal-to-noise, deviations 
from the expected ratio can be attributed to faster or 
slower motion du e to the ecce ntricities of the bodies. 

From Eq 29 in ( |Winn||2010] ), we find that the well de- 
termined ratio is: 



1 + et, sincjb 

1 + Cr sin LOr 



1 



5/3 



(3) 



where the transit parameters that are well determined 
are the periods (P), depths (S), duration (T), and of- 
ten ingress/egress times (r). This ratio is often very 
well determined assuming there is no degeneracy be- 
tween limb darkening coefficients and impact parameters 
(which is not necessarily true for Kepler planet); in the 
case of small eccentricities we expect a good constraint on 
Ch sinwb — Cc sinwc- If this value is large, then it is likely 
that one of the planets is eccentric. If p^. is well known 
from asteroseismology ( Chaplin et al.|2010 Christensen- 



Dalsgaard et al. 2010), then esinw can be inferred for 



each transiting planet (in single or multi-transiting sys- 
tems) . TTVs are also quite sensitive to eccentric ities and 
orientations ( |Nesvorny|2009[ [Payne et aL]|2010[ ) 

Furthermore, if there are multiple large close-in plan- 
ets, then it is reasonable to assume that these planets 
must have migrated inward, perhaps together. Theo- 
ries of multiple-planet migration would p redict that these 



are probably in or near resonances (e.g., Terquem & Pa- 



paloizou 2007), and TTVs are much more sensitive to 



near-resonant planets than ra dial velocities (Agol et al. 
2005 Holman fc Murray|2005 ). Multi-transiting systems 



will be the best for understanding these TTV signals 
without degeneracies, which are present even in well 



chara cterized resonant systems (iMeschiari & Laughlin 
20101. The presence of additional planets with known 



periods also limits eccentricities (past and present) due 
to stability constraints. 

4. A FULL NUMERICAL-PHOTOMETRIC MODEL 

To better understand the value of multi-transiting sys- 
tems, we have developed a model that combines an n- 
body integration with calculation of the photometric 
diminution from each of the planets. This full numerical- 
photometric model can be applied directly in the analysis 
of light curves with interacting planets, including multi- 
transiting systems like those discovered by SIO. While no 
such analysis is performed in this work, we will comment 
on the value of this model in fitting actual data. 

We use the coordinate system suggested by Fab- 
rycky 2009. The orbits are described using astrocen- 
tric equinoctial elements, which avoid the singularities 
of Keplerian elements. The equinoctial elements are 



(Arsenault et al. 1970) a, h — esinnj, k — ecosw, 
p = tan(i/2) sinfl, q = tan(j/2) cosfi, and A™' — M + w, 
where a is the semi-major axis (the same in Keplerian 
and equinoctial elements), e is the eccentricity, zu = u+Q 



is the longitude of pericenter, i.e. the sum of the argu- 
ment of pericenter (w) and the longitude of the ascending 
node (ri), i is the inclination (with respect to the plane of 
the sky), and A™' is the mean longitude, i.e. the sum of 
the mean anomaly (M) and the longitude of pericenter. 
Note that h and k are sometimes defined in the opposite 
way and that A™' is not related to the angle A associated 
with the Rossiter-McLaughlin effect (discussed in [5.2). 



The longitude of the ascending node, f2, in binary stel- 
lar systems is defined as the position angle of the point 
where the companion pierces the plane of the sky com- 
ing toward the obser ver and is me asured with respect to 
astronomical North (Batten 1973). This angle can only 



be determined if the system is somehow resolved on the 
sky, e.g., via astrometry, interferometry, polarimetry, or 
direct imaging. Generally speaking, this angle will not be 
known, and without loss of generality we set Q of the in- 
nermost planet to be equal to 0° . The difference between 
the nodal longitudes of the planets is detectable, so il for 
the remaining planets remains a fitted parameter. 

It is possible that the innermost planets are bright 
enough to be detected in secondary eclipse or in the full 
orbital phase curve. For simplicity, we ignore the phase 
curve since it has more free parameters, but in pr actice 
it is not difficult to include in a parametrized sense ( Kip 



ping fc Bakos 2010). We also do not account for the 
ellipsoidal oscillations of the parent star, which may be 
detectable and can also give c onstraints on the mass ra- 
tios (e.g., Welsh et al. 2010). We choose to keep the 
possibility of fitting the secondary eclipse, since these 
eclipses have significant power to constrain the orbital 
configuration, most notably the eccentricity (Ragozzine 
|fc Wolf 2009) . Of course, they also contain some m- 
formation about the atmospheric reflection and thermal 
properties of the planet, but this is not our focus. Thus, 
our model includes a parameter d'"^'^ which defines the 
depth of the secondary eclipse; if optical and infrared 
secondaries are both observable, these should each have 
their own d'^'^'^ parameters. We note that with current 
techniques, secondary eclipses can only be detected for 
planets with periods of <10 days; planets with larger or- 
bital periods have such weak flux that even Kepler or 
Spitzer photometry cannot detect them. 

The total set of parameters needed to gen- 
erate photometry of a multi-transiting sys- 
tem of two planets is GM*, i?*, ci, C2, C3, C4, 
ai /i?* , /ii , fci , pi , gi , Al , Ml /M* , i?i /i?* , dl"" , 
a2/R*,h2, fc2,P2, 92, A2, M2/M*, i?2/i?*, and, d2'"', where 
all the orbital elements are the osculating astrocentric 
equinoctial elements at a specified epoch (ideally in the 
middle of the dataset) and where generally fii is assumed 
to be fixed at 0° so that pi = tan(ii/2) sin f^i = is 
held fixed. The Ci values are limb-darkenin g coefficients, 
following the notation of Mandel & Agol (|2002|), which 
may not all be used depending on the quality of the 
data. Each additional planet adds 6 orbital parameters, 
2 physical parameters (mass and radius ratio), and a 
parameter for the secondary eclipse depth. 

One of the first steps of any data analysis is to first 
get a good guess of the initial values of all the rele- 
vant parameters. In this sense, multi-transiting sys- 
tems have the nice property that fitting each planet in- 
dividually to the data gives excellent first estimates of 
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nearly all of the parameters described above, especially 
if photometry and radial velocities are both used; this 
is just the normal process used to determine the char- 
acteristics of individual transiting planets. Without ra- 
dial velocities, there is no constraint on the planetary 
masses, but a reasonable range of densities can be em- 
ployed based on the oretical mass-radius relations (e.g., 
Fortney et al. 2007). The main parameter that is not 
well described in this first step is ^^2, the mutual nodal 
angle (as J7 = 0°), which is unconstrained without in- 
cluding the interactions between the two planets. Recall 
that the mutual nodal angle must be determined to cal- 
culate the true mutual inclination between the planets, 
cos (I>i2 = cos ii cos i2 -I- sinzi sin 12 cos(r2i — fl2)- 

Combining the data from both planets to get a global 
fit can be done in two different steps. First, the full pho- 
tometric light curve can be analyzed returning metadata 
like central transit times and transit durations, which 
are then fitted with a fast few-body integrator to re- 
fine the relevant parameters, particularly the best-fit 
masses and Af2. Transit times are most easily deter- 
mined using a numerical integration of the three-body 
problem, though non-resonant systems can be analyzed 
much faster usin g the perturbation method of |Nesvorny| 
& Beauge (2010). After estimates for these parameters 
are narrowed down, the more computationally costly full 
numerical-photometric model can fit the data directly. 
Even when data products like transit times are fit first, 
the full potential of the data can only be realized with a 
numerical model, and it should ideally be the final step in 
characterizing a multi-transiting system. The drawback 
to the full model is that it is much more computation- 
ally time-consuming, especially considering the need to 
explore a highly non-linear ~24-dimensional parameter 
space. 

Even so, a previously unconsidered aspect of multi- 
transiting systems may require a full numerical- 
photometric model to correctly retrieve the desired pa- 
rameters in certain systems. Recall that transit timing 
variations are the result of the aperiodic motion of the 
relative star-planet position in the orbital plane. When 
the relative star-planet position moves out of the orbit 
plane, transit duration variations results from the change 
in impact parameter. Similarly, the relative star-planet 
velocity can change, also causing changes to the shape of 
the transit, though this effect is generally weaker. Each 
of these photometric deviations are approximations to 
the actual case by assuming that the star experiences 
no acceleration during the transit. In certain transiting 
systems, this may not be a valid assumption. An ex- 
treme example is given by the HAT-P-13 system. The 
expected du ration of the puta tive transit of HAT-P-13c 
is 1.28 days ( |Bakos et al.|2009[ ). This is 44% of the orbital 
period of planet b, which is causing the star to executive 
nearly half of its small orbit around the star-b barycenter 
during the transit of c. This affects the relative position 
of the star and planet c at the 1% level, which is creates a 
potentially detectable signature in Kepler-quality data. 
Clearly, the approximation that there is no acceleration 
during the transit is lost and in such a case, the tran- 
sit timing and shaping anomalies are not entirely well- 
defined. Only a full numerical integration of the motion, 
coupled directly to the photometric model, will be capa- 
ble of a precise description of the data. We note that 



this new aspect of multi-transiting systems gives some 
additional constraint, albeit weak, on the orbit of one 
planet from its induced stellar motion during the transit 
of other planet that may not be entirely captured by a 
TTV/TDV analysis. 

To account for these and other small effects (such as the 
slightly asymmetric nature of eccentric transits and light- 
travel time effects), a full numerical-photometric model is 
employed. The model can also easily calculate stellar ra- 
dial velocities and eventually the Rossiter-McLaughlin ef- 
fect. In our implementation of such a model, the equinoc- 
tial elements are converted to Cartesian coordinates and 
passed to a Fortran 90 n-body integrator which inte- 
grates the few-body system. This integrator can also be 
readily expanded to include non-Keplerian effects from 
General Relativity, and the rota tional and tidal bulge s 
of the star and planet (following 'Mardlin g fc Lin||2002 ), 
though these are not included in the present analysis. 

There can be a difference between the positions of the 
star and planets in the timeframe of the system and their 
positions as observed by the Kepler satellite, due to vari- 
ations in the light travel time. First, we assume that 
the motion of the Earth or telescope has been correctly 
accounted for, i.e. all measurements are referred to the 
barycenter of the solar system (BJD). Any remaining off- 
set is due to the motion of the star around the barycen- 
ter of the system and, perhaps, motion of the system 
barycenter relative to the solar system. This latter ef- 
fect can have an observable consequence in certain cases 
( |Scharf||2007| lRafikov||2009[ ) , but is not considered here. 

W e found tha t accountmg for the light-travel time ef- 
fect ( |Loeb|2005D was best accomplished by correcting the 
positions of the star and each planet to the plane of the 
sky at the system barycenter at a consistent time. These 
positions are then used to calculate photometry. 

Instantaneous stellar velocities (for comparison to ra- 
dial velocity data) can be calculated from the positions. 
For all planets, i, all the times at which the on-the- 



yf is less than the separation 



sky distance di 

needed for transits (i?j + R^,), the photometric di mming 
due to the transit is determined from the codes of lMan-l 
del & Agol (2002[), using an appropriate limb darkening 



law and coetticients. 

The quantity of Kepler photometry is in a new regime 
of data fitting that suggests a slightly alternative method 
for calculati ng photometric light curves. The codes of 



Mandel fc Agol| (|2002j) are fast, but not instantaneous, 
so using them less often can improve computation time. 
We find that it is much more efficient to calculate the 
intensity as a function of planet position (with given 
limb- darkening coefficients and radius ratio) with a sam- 



pling of about 200 positions from to 



R».-\- Ftp 



Fast cu- 



bic interpolation is used to calculate the intensity at all 
the thousands (or more) of planet positions needed for 
the observations. One could imagine precomputing all 
the intensity calculations and then using a large multi- 
dimensional look-up table to calculate photometry, but 
we concluded that the interpolation method we employed 
is sufficient. 

The model flux is calculated every ^100-1000 seconds 
(the time step of the integrator) and then an integration 
of the cubic interpolate is used to calculate the intensity 
given the beginning and end times of each observation 
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(Kipping 20101. Note that this correctly accounts for 
tlux changes that may happen within a single integration 
(Kepler long cadence observations which have 30 minute 
long integrations) and can also account for any partic- 
ular duty cycle (including aperiodic observations, closed 
shutter time, read-out time, etc.). We have verified that 
the timestep size and interpolation step sizes are small 
enough to not affect the light curve at the 1 part per mil- 
lion level or (usually) less in a variety of circumstances. 
The photometry and the interpolation is only performed 
near transits and eclipses to speed up the calculations. 

The flux diminutions due to all the planets is added 
together in order to calculate the model brigh tnes s of 
the system at the observed time. (See Section 5.1 be- 
low for a discussion of what happens when the planets 
overlap each other during transit.) We are also assum- 
ing that slow variability due to the star or other data 
reduction trends on timescales longer than the transits 
have been successfully removed. Our method, like other 
methods, also assumes there are no significant changes 
in brightness during a single photometry calculation of 
^100 seconds. This is true for transit and eclipse light 
curves, but can be more complicated for certain kinds of 
fast mutual events. For this reason, mutual events are 
not typically included in our photometric calculations. 

When fitting data, the residuals of the model from the 
data are computed, squared, and summed to calculate 
in the usual fashion. Data far outside of transits or 
eclipses is not used in the calculation of x^- For ~100- 
day long simulated Kepler photometry at short cadence 
with a precise tolerance, a single computation takes 
approximately 1 second. To estimate the most likely val- 
ues for our parameters now requires a minimization of 

■ This highly non-linear multi-dimensional parameter 
estimation is difficult, but is assisted significantly by the 
fact that the independent Keplerian fits to the transit 
photometry and RV data give an excellent first approxi- 
mation of most of the fitted parameters. 

An example of some of the simulated photometry from 
this model is given in Figure [3j Using this model to 
fit simulated photometry, we have drawn the following 
conclusions. 

• The global parameter space is covered with many 
disconnected local minima. We have deve loped 
an MC MC- like c o de (b ased on concepts in Fbrd 
j2005| ) and |Ford| ( |2006[ )) that is coupled witlTa 
Levenburg-Marquardt minimization that is run ev- 
ery ^50 links (but this is not used to update the 
MCMC chain as this would violate the Metropolis- 
Hasting algorithm). This combination of a local 
and global minimizer to explore parameter space 
and find the best minima is also employed by 
Meschiari & Laughlin (20101. When started near 
the global minimum, this code is successful at find- 
ing good solutions. For real data, it will be diffi- 
cult, if not impossible, to be sure that the code has 
settled in the best global minimum, although re- 
duced values near 1 with appropriate errors are 
a good sign. Estimating reasonable errors on the 
system parameters is also very difficult, especially 
with limited computational time. 

• The detectability and precision with which param- 
eters are estimated is in accordance with the ex- 




400 600 
Days 



1000 



Fig. 3. — Simulated Multi- Transiting System showing TTVs. 
Using the numerical-photometric model described in the text, the 
light curve shown represents two hot Jupiters in the 2:1 resonance 
with periods of approximately 3 and 6 days and near central impact 
parameters. This model includes self-consistently the three-body 
interactions, light travel-time corrections, and the acceleration of 
the star from one planet during the transit of another planet, which 
is not entirely captured in models that use only TTVs and TDVs 
(see text). The light curve covers 1000 days of simulated Kepler 
photometry at 30 minute ("long") cadence with no noise added. 
The two planets were given different depths; the outer planet (with 
less frequent transits) is clearly larger. Occasionally both planets 
transits the star at the same time, these double transits have a 
combined depth grea ter t han either planet individually; mutual 
events (described in i ]5.1| and Figure H are not included in this 
model. The first two double transits nappen coincidentally and 
such events would repeat about every year in the purely Keple- 
rian case. The train of double transits around day 750 is caused 
by planet-planet perturbations and changes in the nominal orbital 
frequencies and demonstrates that transit timing variations are in- 
cluded in the model. 



pected results. For examples, in systems where 
dynamical interactions are important and where 
TTVs would be observed, we recover the correct 
value of the mass ratios of both planets and we 
determine the semi-major axes very precisely (as 
seen also in Nesvorny & Morbidelli 2009). If com- 
bined with precise radial velocities, as in Meschiari 
& Laughlin 2010, one could measure dynamical 
masses of each of the objects. 



• §3.2| discussed the possibility of measuring dynami- 
cal masses from photometry alone by observing the 
light-travel time effect. Even with 1000-days of Ke- 
pler short-cadence photometry, this will not be pos- 
sible in most systems. In a system with two tran- 
siting interacting hot Jupiters, the true mass of the 
system was only constrained to within a factor of 
~2 with 100-days of data. This implies that (typ- 
ical) ground-based follow-up observations will be 
important for understanding multi-transiting sys- 
tems. 

• Since the two planets generally cross the star along 
two different chords, one might imagine that multi- 
transiting systems are efficient at breaking limb- 
darkening parameters. Indeed, this is a major 
motivation for creating a single model that in- 
corporates data for both planets simultaneously. 
Even so, we found that at the precision expected 
from Kepler even on relatively faint stars (1000 
ppm/min, V~14), the limb darkening parameters 
are already measured by the largest planet and the 
addition of another planet does not substantially 
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improve the measurement of limb darkening. How- 
ever, this conclusion was not tested in detail. 

5. OTHER TECHNIQUES FOR STUDYING 
MULTI-TRANSITING SYSTEMS 

5.1. Mutual Events 

Systems with more than one planet have the unique 
possibility of photometric mutual events, where the ob- 
jects shadow or occult one another as seen from Earth. 
The primary transit and secondary eclipse are well- 
known versions of these mutual events, since these "star- 
planet" mutual events are the only phenomena relevant 
for single planets. When additional bodies are added, 
however, other mutual event phenomena are possible. In 
multi-transiting systems, the likelihood of observing and 
correctly interpreting these mutual events is much higher 
than in systems where one or no planets transit. Fur- 
thermore, if they can be observed, mutual events pro- 
vide unique constraints on the orbital parameters of the 
planets involved and are worth investigating in detail. 

As in j ]2.1| above, we begin by being clear with termi- 
nology. In a multi-transiting system, when two planets 
transit at the s ame time , this is called a double transit 
(following, e.g., Holman 20101. If a mutual event be- 
tween the planets occurs during the transit, we call this 
an overlapping double transit, which leads to a bright- 
ening in the light curve as described below. When two 
planets cross in front of one another outside of transit, 
we call this a planet-planet occultation. Planet-planet 
shadowing occurs when an inner planet passes between 
an outer planet and the star. Finally, following termi- 
nology from the solar system, a mutual event between 
three objects in an extra-solar system could be called a 
syzygy. Other than the overlapping double transit case 
(a star-planet-planet syzygy), we do not consider other 
exo-syzygies here. The systems of SIO do not show any 
double transits in the 43 days of data shown. 

Mutual events are very powerful in constraining the 
orbital and physical properties of the planets. When a 
mutual event occurs, it is like a precise differential astro- 
metric measurement and one can know almost exactly 
where both of the participating planets are at a specific 
time. In particular, we find that mutual events are ex- 
cellent at measuring the mutual nodal angle, and hence 
the true mutual inclination, between the participating 
planets (see discussion below). 

In systems were many mutual events can be observed, 
since these events involve intersections of comparatively 
small planets, they result in precise requirements that 
must be correctly reproduced in a viable model. One 
can imagine tracking "mutual event timing (and dura- 
tion) variations" which are even more sensitive to small 
orbital perturbations than transit timing variations and 
consequently provide even stronger constraints on poten- 
tial perturbations. 

We have considered all the possible mutual events be- 
tween a star and two planets and here describe the pos- 
sible mutual event phenomena in order of detectabil- 
ity: overlapping double transits, planet-planet occulta- 
tions, and planet-planet shadowings. Note that these 
planet-planet mutual events are intrinsically more rare 
and shorter-duration than star-planet transits since plan- 
ets are much smaller than stars (and they are not in orbit 
around one another). 



5.1.1. Overlapptng Double Transit 

If two planets overlap while they are in transit, then 
the photometric diminution is no longer the sum of two 
individual diminutions due to the two planets. The shape 
of the covered part of the star deviates significantly from 
two circles and standard light curve calculations can no 
longer be used. Even when both planets are known to 
transit, the probability of such an event depends on the 
mutual inclination and node of the planets and these 
events can occur often or never, even if (non-overlapping) 
double transits happen regularly. Even if the planes of 
the planets are properly aligned, the phases of the plan- 
ets often make this a rare occurrence in time (assuming 
resonances are not in play). 

Even so, when an overlapping double transit does oc- 
cur, it can have easily detectable photometric signal. As 
mentioned above, early observations by HST of TrES-1 
saw a bump in the light curve that was tho ught to per- 
haps be due to a overlappi ng double transit ([Rabus et al. 
2009 1 , but later analysis ( Dittmann et al. 2009 ) showed 
that a starspot cross ing was a much more likely conclu- 
sion (see also ^ 5.1.4 below). When the orbits of both 
planets are known well, and given the exquisite knowl- 
edge of the periods and phases of high signal-to-noise 
transiting planets, there will be much less confusion on 
this point. Overlapping double transits are also the eas- 
iest mutual event to identify: they can only occur dur- 
ing a double transit and are characterized (in the best 
cases) by an unmistakable brightening (Figure |4]) . Plan- 
ets crossing over starspots can also cause minor bright- 
enings, but the overlapping double transit usually has 
enough difference in amplitude and duration that we find 
it unlikely that a starspot crossing will be confused with 
a mutual event, especially if the event is not grazing. 

If these events happened over a source of uniform 
brightness, the expected lightcurve could be generated 
analytically. As the star is limb-darkened with both 
planets moving across it and overlapping, we find that 
a numerical technique for modeling the light curve is re- 
quired. Since these events are generally rare, this is not 
computationally costly overall and such a model can be 
folded into the numerical-photometric model described 
above. To correct the light curve created from the sum 
of the two individual light curves requires adding back in 
the flux that was doubly subtracted from the region of 
overlap between the two planets. We have generated an 
efficient method for calculating this light curve (similar to 
Barnes & Fortney 2003) and tested it with Monte Carlo 
integrations (similar to Carter & Winn 2009). Resulting 
light curves of overlapping double transits are given in 
Figure Hp. 

We have modeled the case of an overlapping double 
transit for two Neptune-size planets in circular orbits 
around a relatively faint Kepler star {V ~ 14) and find 
that light curves for mutual inclinations that differ by 
only a few degrees are distinguishable. Indeed, since 
both planets are transiting the star, all the orbital param- 
eters are relatively well known and the only parameter 
with significant freedom is f2i — 112, the mutual nodal an- 
gle; since this parameter strongly affects the timing and 

^ An animation of an overlapping dou- 
ble transit is available for download at 
http: //www. cf a. harvard, edu/^dragozzi/meanim. gif . 



Multiple Transiting Planets 



13 





1 


000 ; 


X 





999 r 








ive F 





998 - 


D 





997 r 


q: 





996 [ 







995 L 




200 400 600 800 1000 
Time (minutes) 

Fig. 4. — Simulated Overlapping Double Transit. If two tran- 
siting planets were to overlap while crossing the stellar disk, a mu- 
tual event we call an overlapping double transit, they would cover 
a smaller area of the star and a brightening would result. This 
is shown in the above figure with a simulated mutual event, the 
observed photometry is given by the solid line and would have a 
distinct "W" shape. The dashed line shows what the transit of the 
outer planet would look like alone and the dotted line (mostly con- 
current with the solid line) shows what the double transit would 
look like without a mutual event. The difference is a significant 
brightening whose signal is almost as big as the inner planet tran- 
sit, since the overlap is nearly complete. Inset is a picture showing 
the two orbital paths (straight lines), their mutual inclination of 5 
degrees, and three positions of each planet to show which orbital 
path corresponds to which object (they are moving from right to 
left). The filled figure is the coverage of the two planets at the 
moment of closest approach. Although these are only Neptune size 
planets, a mutual inclination different by only a few degrees would 
have a detectably different signature by Kepler even around a faint 
star. Mutual events are unique and powerful ways of constraining 
the true mutual inclination, as described in the text. An anima- 
tion of an overlapping double transit is available for download at 
http: //www. of a . harvard . edu/^dragozzi/meanim.gif . 



shape of mutual events, it is not surprising that these 
events are very good at determining f2i — ^2 without 
degeneracy. Other parameters, like the eccentricity, can 
also be constrained as they also affect both the exact po- 
sitions and exact speeds, which are measured precisely 
by mutual events. 

A simple geometric argument can explain why this is 
the case. Consider the detection of an overlapping double 
transit, where fli is set to 0° without loss of generality. 
This means that the inner planet traverses an essentially 
horizontal line across the star. The impact parameter is 
known and the time of the mutual event gives the ex- 
act orbital phase: from the remainder of the transit light 
curve, this phase is directly related to the planet's po- 
sition relative to the center of the star and the exact 
position of the inner planet is known (modulo unimpor- 
tant rotation in the plane of the sky). The second planet 
is located in the same position and, from its impact pa- 
rameter and phase, is a specific distance from the limb 
of the star, i.e., the locus of possible limb crossings is a 
circle of fixed radius (since its orbit can be oriented in 
any direction, not just horizontally). This circle inter- 
sects the actual limb of the star in one or two locations, 
immediately limiting the possible trajectories to fewer 
than two possibilities, which translates into a two-way 
degeneracy in the true mutual inclination. Including in- 
formation of the shape of the mutual event light curve is 
generally sufficient to break this degeneracy; an extreme 
example of this occurs when one planet is prograde and 



the other is retrograde, where the mutual event will be 
much shorter than in a system where both planets orbit 
in the same direction, due to the difference in apparent 
relative velocity. Observable out-of-transit planet-planet 
occultations have a similar effect of measuring precise 
orbital parameters, including mutual inclination. 

5.1.2. Planet-Planet Occultation 

Even in multiple-planet systems that are not transit- 
ing, it is possible for the planets to c ross one another as 
seen from the perspective of Earth. Cabrera & Schnei- 
der ( 2007 1 consider such events (as well as planet-planet 
shadowmgs discussed b elow) for p lanets and large moons 
or bina ry pla nets as do Sato fc Asada^ (2009) and Sato Sz 
Asada (2010). These authors discuss such niutual events 
at length tor both optical and thermal measurements, but 
do not discuss mutual events between unbound objects 
(e.g., two objects in independent astrocentric orbits). 

The frequency of planet-planet occultations type of 
event depends significantly on the geometry of the sys- 
tem. If the planets are coplanar with each other and at 
a very low relative inclination to the Earth (e.g. nearly 
central impact parameter) , then occultations could occur 
very frequently. The probability of such an alignment is 
about — , which is <0.5% even in the most ideal case. 
On the other extreme, planets in very different planes 
can have rare and aperiodic events. In the case of cir- 
cular orbits, the apparent orbit on-the-sky is an ellipse 
with major axis a and minor axis acosi, with the minor 
axis of the ellipse oriented at an angle of $7 east of astro- 
nomical North. Mutual events are only possible if these 
apparent orbits have at least one point where the mini- 
mum distance between them is the sum of the planetary 
radii. This will be most likely for systems with different 
values of $7, i.e., those with at least some mutual incli- 
nation. Though non-transiting planets can participate in 
planet-planet occultations, predicting and understanding 
these events is far easier in multi-transiting systems. 

Planet-planet occultations are very weak, since it is 
the reflected/emitted photons of the blocked planet that 
are being lost. They are completely unobservable in sys- 
tems where the light from the occulted planet cannot 
be detected, which limits its applicability to occulta- 
tions of (large) planets with periods <10 days (where 
reflected light or thermal emission is detectable at some 
wavelengths with current techniques). None of the sys- 
tems in SIO are likely to show planet-planet occultations. 
The duration of these events is also much shorter than 
the secondary eclipse, the usual source of occulting the 
planet, since the physical sizes of the objects are so small, 
though the velocity of the planets on the sky may in un- 
usual cases extend the length of such events. Even in 
the most favorable cases, durations are no longer than 
^•1 hour. Even with fantastic photometry from Kepler, 
Spitzer, or Hubble space telescopes, individually detect- 
ing these events is difficult even when the outer planet 
occults a hot and/or reflective inner planet. However, if 
such events could be detected (or ruled out) with rea- 
sonable signal-to-noise, they would place very tight con- 
straints on the orientation of the planetary system, since 
there would be information on the relative astrometric 
position along a line of sight different than that con- 
necting the Earth the star (see the geometric argument 
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above) . 

Furthermore, in these events, high signal-to-noise ob- 
servations can even characterize the reflected/emitted 
hght spatially on the occulted planet, possibly leading 
to surface maps that are not (mostly) degenerate in lati- 
tude like those pr oduced by phase curv e s and secondary 
eclips e mapping (Knutson et al. 2007 Rauscher et al 



2007 1 . We n ote that this is al so true m overlapping dou- 
ble transits (Kipping & Tinetti 2009), though the signal 



due to the planetary occultation in that case is much 
more difRcult to pick out. 

It might be possible, in well-characterized systems, to 
predict such events and observe them in the infrared, 
where the signal is likely to be higher. The James Webb 
Space Telescope may be able to detect some of these short 
but highly interesting events. 

5.1.3. Planet-Planet Shadowing 

The smallest effect we consider is planet-planet shad- 
owing. As planets orbit, the shadow of planet 1 can cross 
planet 2, reducing the amount of light it reflects; from 
the perspective of outer planet, it is observing a transit 
of inner planet. If the two planets are distant from one 
another, the reflected light from shadowed planet will 
drop by the depth of the transit, i.e. ^1% if the in- 
ner planet is Jupiter-like. Since the reflected light signal 
is so faint anyway (10~^ at best), this becomes a part 
per million signal that is very small. You do get a geo- 
metric enhancement of about 1 H — — , but even in the 

02— ai ' 

closest orbits, this only amounts to an enhancement of 
a factor of a few. This effect is maximized for "double 
planets" as has been discussed in this conte xt by [Cabrera 



& Schneider! ( |2007[ ), |Sato &: Asada| ( |2009[ ), and [Sato fc 
Xsada (20OT. 

On the other hand, planet-planet shadowings are likely 
to be present in systems with low mutual inclinations, 
they repeat at the potentially-frequent synodic period of 
the planets, and the duration of such events can be rel- 
atively long. So, it might not be impossible to detect 
this effect and it could even give constraints on the or- 
bital properties of the planets. However, since the signal 
is so small, observing planet-planet shadowing may not 
increase our knowledge of the system if both planets are 
transiting anyway, though it may be a systematic source 
of small residuals. Modeling this effect requires detailed 
knowledge of the reflected light curve of the shadowed 
outer planet and would require correct modeling of other 
effects (such as variability due to starspots or planetary 
weather) . 

5.1.4. Other Effects Related to Mutual Events 

Mutual events may be relatively rare; very rough es- 
timates indicate that Kepler should (serendipitously) 
observe a few overlapping double-transits, though this 
likelihood depends strongly on the actual prevalence of 
multi-transiting systems. It is more common to expect 
transits of planets across starspots (which have already 
been observed more than once in known transiting sys- 
tems). In addition, Kepler and CoRoT provide long- 
term light-curve fluctuations that can be used in con- 
junction with starspot t ransit s, to construct spot mod- 



use for determining orbital parameters may be limited. 
Besides an under-constrained spot model, trying to in- 
fer the mutual inclination from starspot crossings also 
requires some knowledge of the stellar spin axis and dif- 
ferential rotation, though perhaps these limitations can 
be overcome in some systems. 

A few very minor effects are worth noting here for com- 
pleteness. If a close-in planet reflects light back onto the 
star or otherwise creates a subplanetary hotspot, this 
can have an effect on the transit light curve of other 
planets. Similarly ellipsoidal oscillations (and the associ- 
ated gravity darkening) due to one planet can change the 
transit properties of another. Very close-in planets with 
escaping atmospheres could have different mutual event 
properties in different filters, which might help probe the 
shape of the escaping atmosphere. 

5.2. Radial Velocities and Multi-Rossiter-McLaughlin 

Though a surprising amount of information about 
multi-transiting systems can be inferred from photome- 
try alone, radial velocity measurements will improve and 
confirm other aspects of these systems. As mentioned 
in |3.2[ the combination of TTVs from photometry and 
RVs gives dynamical masses and absolute radii, which 
is very exciting. RV observations are also important to 
search for additional non-transiting planets that do not 
create observable TTVs. 

Radial velocity observations during transits have par- 
ticular value: as a transiting planet crosses the parent 
star, an anomalous radial velocity signal is generated 
due to the occultation of portions of the rotationally 
Doppler shifted star. This effect is known as the Rossiter- 
McLaughlin effect and is sensitive to the projected an- 
gle between the stellar spin axis and the orbital pole of 
the planet, called A (equivalent to — /3). This effect has 
been measured for several transiting planets (Triaud et 
al. 2010), with some planets nearly co-planar with their 
corresponding stellar equators and many planets showing 
p olar or possibly r etrograde orbital paths. 



[Fabrycky ( |2009 1 discussed the possibility of using the 
RJVL effect to h elp constrain the p lanetary mutual incli- 



els, as in Huber et al. 



tems, this ability is o 



(|2010|). In multi-transiting sys- 
y pronounced, though its 



nation ((/)) and Simon et al. ( |2010 ) consider the value of 
RM on understanding the orbit of transiting exomoons. 
Here we expand on this discussion to be more specific 
about the importance of "multi-RM" observations. 

Note that the orbital normal to each planet and the 
spin vector of the star all can be described by a single unit 
vector or by 3 Euler angles, assuming that the starting 
vector is pointing in the z-direction. Using w*, i*, and il* 
as Euler angles in the standard way for three-dimensional 
orbit description proscribe the stellar spin vector as well. 
(For rotational axes, w* is not important, unless you're 
interested in the rotational longitude at epoch.) 

Using this notation, the true angle between the star's 
spin axis and a planet's orbital axis, ■0, can be written 
as 

cos "0 = cos cos ip + sin sin ip cos(il* — ilp) (4) 

and comparing this to the definition of A from Fabrycky] 
|fc Winn (2009) , we find that A = fi* — fJp, i.e., the differ- 
ence m nodal angles. Whereas in the case of transiting 
planet orbits, where both inclinations are known, but the 
relative nodal angle is not, the Rossiter-McLaughlin ef- 
fect is sensitive to the relative on-the-sky alignment angle 
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between the star and planet, but the stellar inclination 
to the line-of-sight, is not determined. 

This difference can be exploited in multi-transiting sys- 
tems and it is trivial to show that fli — — \Xi ± A2I. 
This means that if two transiting planets have measured 
Rossiter-McLaughlin angles, the true mutual inclination, 
4>, is determined to a two-way degeneracy. If one of Ai 
or A2 is small (i.e. aligned with the star) then the de- 
generacy is weak and mutual inclination can be deter- 
mined well, assuming that the errors on the values of A 
are small (often the errors are substantial, especially on 
fainter stars). In a triple-transiting system, measurement 
of all three values of A lift the degeneracy completely, and 
all the mutual inclinations are directly measured. 

This is a promising result, since the other major meth- 
ods for measuring true mutual inclination - TTVs and 
mutual events - will not be possible for all systems 
and/or may require many years of observations. Multi- 
RM is typically possible in all multi-transiting systems 
and is limited only by the brightness of the star, stellar 
type, stellar rotational velocity, the sizes of the planets, 
and other observational constraints. This suggests that 
a campaign of multi-RM measurements can be used to 
probe the mutual inclination distribution of planets and 
provide imp ortant insights into their formation and evo- 
lution (^3.3 1. Such a campaign would also, clearly, probe 
projected star-planet alignment for each of these planets, 
which is also relevant to their formation and evolution 
and should provide useful clues to understanding various 
properties of star-planet alignment in general (Triaud et 
al. 2010; Winn et al. 2010, submitted). 

Finally, since two different planets cut different chords 
across the star, it may be possible to break some of 
the degeneracies normally associated with the RM ef- 
fect, though the effects of differential rotation should be 
considered in this regard. For example, the (latitudinal) 
motion of star spots depends on the unknown , and the 
timing o f spot crossings can pote ntially constrain this pa- 
rameter ( Miller- Ricci et al.|2008| ; the likelihood of useful 
spot crossings seems much nigher in multi-transiting sys- 
tems and Kepler provides long-term light-curve fluctua- 
tions that can be used to constrain spot models (S 5.1.4). 



6. CONCLUSIONS 

The announcement of candidate multi-transiting sys- 
tems by SIO ushers in a new era for transiting exo- 
planets. Besides great potential for informing theories 
of planetary formation, planetary interiors, orbital mi- 
gration, large-scale scattering, and other effects, multi- 
transiting systems are among the best for understand- 
ing the question of planetary habitability. The radii of 
potentially habitable super-Earths are measured photo- 
metrically, and their masses can sometimes also be mea- 
sured photometrically by TTVs induced on other tran- 
siting planets, which is especially valua ble in the cas e 
where RV measurements are not feasible ( Holman|2010 ) . 
Only with good estimates for masses and radii can hab- 
itability be seriously assessed. Furthermore, the orbital 
architecture of the entire system is much better under- 
stood if more than one planet is transiting the star, which 
is important for such questions as long term stability and 
potential climate cycles. 

We have attempted a systematic investigation of the 
value of systems with multiple transiting planets. Many 



of the basi c ideas have been presented by other authors 
(especially 'Schneide r|200H |Fabrycky|2009{ |Holman|2010 
|Steffen et al. 2010) , but we have expanded their results 
in a way that will be helpful for maximizing the value 
of the observations. This has resulted in the following 
conclusions. 

1. Based on geometric constraints and RV observa- 
tions, multi-transiting systems are likely to be com- 
posed of close-in multiple Neptune systems with 
low mutual inclinations. While the frequency of 
these systems remains to be seen, the identifica- 
tion of at least five good multi-transiting candidate 
systems by SIO using only 43 days of data implies 
that Kepler will discover several multi-transiting 
systems. 

2. Transiting planets in general are useful for obtain- 
ing important physical parameters about distant 
exoplanets, including mass, radius, and density. 
Multi-transiting systems do this efficiently for more 
than one planet in the same system. Comparative 
planetology between these planets is more powerful 
than between planets in disparate systems. 

3. Multiple planet systems in general are useful for 
understanding the formation and evolution of plan- 
etary systems. Multi-transiting systems are multi- 
ple systems that will be very well characterized, 
with true masses and other orbital parameters 
(such as resonance occupation) precisely measured, 
along with a well-determined orbital architecture 
(although potentially missing some non-transiting 
planets that are too small to generate observable 
TTVs or RVs). 

4. Multi-transiting systems can be mimicked by false 
positives, and each candidate planet should be con- 
firmed separately, though it is increasingly diffi- 
cult to develop realistic astrophysical models for 
systems with many planet candidates. Consistent 
inferred stellar properties (i.e., density and limb 
darkening) are important to check, but are not 
generally sufficient to prove the planetary hypothe- 
sis. Once the candidates can be shown to orbit the 
same star, (e.g., through the identification of robust 
and self-consistent TTV signals), it is generally 
straightforward to prove that their masses must be 
small (e.g., through TTVs (or lack thereof), orbital 
stability, or other theoretical grounds). 

5. Multi-transiting systems that show TTVs are much 
easier to interpret since the size, period, and phase 
of the perturbing planet are known. Combination 
of TTVs with RV measurements can allow for the 
determination of dynamical masses and abs olute 
radii for ea ch of the objects in the system (Agol 



et al. 2005). In the best cases, dynamical mass 
limits can be placed on these systems through pho- 
tometry alone, if li ght-travel time or pho tometric 



Doppler boosting (Loeb & Gaudi 2003) are de- 
tected, but basic ground- based follow-up will al- 
ways be justified. 

6. Combining the inclination distribution of exoplan- 
ets with other properties, such as eccentricity and 
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star-planet alignment distributions, will be ex- 
tremely valuable for interpreting the formation 
and evolution of individual multiple systems and 
planetary systems in general. Although the true 
mutual inclination is not determined directly in 
multi-transiting systems, there are several useful 
techniques that can be applied to these systems 
that can constrain this valuable parameter includ- 
ing TTVs, exoplanet mutual events, and multi- 
ple Rossiter-McLaughlin measurements. Multi- 
RM observations constrain the true inclination 
with comparable errors as the star-planet align- 
ment angle, A, though in cases where two plan- 
ets are mis- aligned there is a two-way degener- 
acy. Multi-transiting systems offer the best hope 
for probing the inclination distribution of plane- 
tary systems. 

7. We have developed a full numerical-photometric 
model that accounts for all the important effects 
of interacting planetary systems and photometric 
diminutions. In such a model, the acceleration of 
the star due to a perturbing planet can be mea- 
sured and characterized; this is not entirely cap- 
tured in models that only fit TTVs and TDVs. The 
model also exploits the fact that both planets tran- 
sit the same star (with the same radius and limb 
darkening), and accurately accounts for light-travel 
time offsets. Due to computational limitations, this 
model is best used as the final step in determining 
orbital and physical parameters in multi-transiting 

Pulsar timing can measure orbital characteristics of "pul- 
sar p lanets" much mo re precisely than photometric measurements 
(e.g. , |Wolszc2an|1994[| . However, the radii of pulsar planets are not 
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8. Mutual events between exoplanets, including over- 
lapping double-transits, planet-planet occultations, 
and planet-planet shadowings, are very valuable, 
especially for measuring the mutual inclination be- 
tween two transiting planets. The strongest signal 
is given by overlapping double-transits, for which 
we have developed an efficient model. Though not 
as precisely constrained as mutual events, starspot 
crossings may also have some value as they can be 
more common. Finally, planet-planet occultations 
may be valuable for probing the surface brightness 
distribution of exoplanets in a unique way. 

We have seen from the study of stars that eclips- 
ing double-lined spectroscopic binaries are the observa- 
tional foundation on which much of stellar astrophysics 
is based. We propose that stars with multiple transiting 
planets will take a similar role in characterizing extra- 
solar planets and shaping our models of the formation 
and evolution of planetary systems. For the reasons out- 
lined in this work, multi-transiting systems are the most 
information-rich planetary systems outside our solar sys- 
tem around main sequence stars^, and should receive spe- 
cial attention. 
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